Spin-orbit interaction and asymmetry effects on Kondo ridges at finite magnetic field 
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We study electron transport through a serial double quantum dot with Rashba spin-orbit inter- 
action (SOI) and Zeeman field of amplitude B in presence of local Coulomb repulsion. The linear 
conductance as a function of a gate voltage V g equally shifting the levels on both dots shows two 
B — Kondo ridges which are robust against SOI as time-reversal symmetry is preserved. Resulting 
from the crossing of a spin-up and a spin-down level at vanishing SOI two additional Kondo plateaus 
appear at finite B. They are not protected by symmetry and rapidly vanish if the SOI is turned 
on. Left-right asymmetric level-lead couplings and detuned on-site energies lead to a simultaneous 
breaking of left-right and bonding-anti-bonding state symmetry. In this case the finite- B Kondo 
ridges in the V g -B plane are bent with respect to the V 9 -axis. For the Kondo ridge to develop 
different level renormalizations must be compensated by adjusting B. 

PACS numbers: 05.60.Gg, 71.10.-w, 73.63.Kv, 76.20.+q 
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I. INTRODUCTION 

In linear response transport through quantum dots the 
spin Kondo effect shows up a as a plateau in the linear 
conductance G when varying the level positions by an 
external gate voltage V g , often referred to as a Kondo 
ridge. 1 -"— The width of the Kondo ridge is set by the local 
Coulomb interaction U on the dot, which, in the Kondo 
regime, exceeds the level-lead hybridization T. The lat- 
ter determines the width of the Lorentzian resonance 
in G(V g ) at U = 0. Breaking the two- fold Kramers- 
degeneracy by a local Zeeman field of amplitude B de- 
stroys the Kondo ridge; along the V^-axis the conduc- 
tance plateau is split up into two Lorentzian resonances. 
In contrast, spin-orbit interaction (SOI), although break- 
ing spin-rotational symmetry by designating a certain 
(spin) direction, does not destroy the Kondo effect 4"— 
In the presence of SOI spin is no longer a good quantum 
number but a Kramers doublet remains as time-reversal 
symmetry is conserved. This leads to a Kondo effect in 
the presence of a local interaction provided the gate volt- 
age is tuned such that the dot is filled (on average) by an 
odd number of electrons (dominant spin fluctuations). 

In multi-level dots increasing B might lead to energet- 
ically degenerate states (level crossings) resulting from 
different orbitals. If one is a spin-up state and one a spin- 
down state and the gate voltage is tuned such that an 
electron fluctuates between these states one might expect 
the emergence of a spin Kondo effect at finite magnetic 
fields* — If the orbital quantum number is conserved 
in the leads in such systems additional orbital Kondo 
effects^ and combinations of spin and orbital Kondo 
effects^ may appear. Here we consider a setup where the 
orbital quantum number does not arise in the leads and 
we thus concentrate on spin Kondo effects. In contrast 
to the standard B = Kondo effect the one appearing at 
finite B is not protected by time-reversal symmetry and 
we show that it is destroyed in the presence of a finite 



SOI. 

We here study a serial double quantum dot, each hav- 
ing a single spin-degenerate level (at vanishing magnetic 
field) described by a tight-binding model with two lat- 
tice sites coupled by electron hopping of amplitude t 
and connected to two semi-infinite noninteracting (Fermi 
liquid) leads via tunnel couplings of strength Tr, and 
T r. The on-site energies of the two levels are given by 
e i/2 = Vg i 6. The Rashba SOI identifies the z-direction 
of the spin space and is modeled as an imaginary elec- 
tron hopping with spin-dependent sign between the two 
lattice sites We here exclusively consider the cou- 

pling of a magnetic field to the spin degree of freedom 
(Zeeman term) and neglect its effect on the orbital mo- 
tion. The relevant component of the Zeeman field per- 
pendicular to the SOI defines (without loss of generality) 
the x-direction. The local Coulomb repulsion is modeled 
as an on-site U as well as a nearest-neighbor U' repul- 
sion and treated within an approximate static functional- 
renormalization group (fRG) approach^ Our model is 
sketched in Fig. [1] 

In the absence of SOI a finite- B Kondo effect will lead 
to a Kondo ridge. This conductance plateau can be de- 
tected if G is computed (or measured) as a function of V g 
and B. We show that the finite- B conductance plateau 
develops on a line parallel to the V^-axis, the only excep- 
tion being the case with broken left-right and bonding- 
anti-bonding state symmetry realized for Tr, =/= Tr and 
(5^0. In this case the finite- B Kondo ridges are bent 
with respect to the V^,-axis; to follow the conductance 
maximum when changing V g one has to adjust B to 
compensate for the asymmetry-induced level renormal- 
ization. As the two symmetries are generically broken in 
experimental systems our results are important for the 
understanding of measurements of the linear conductance 
of multi-level quantum dots as a function of Vq and B. 

This paper is organized as follows. In the next section, 
we introduce our double-dot model. In Sec. IIIII details of 
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FIG. 1. The considered setup consists of a serial double quan- 
tum dot with energies = V g ± S coupled by a hopping 
amplitude t and a SOI of strength a. The levels are split by 
an external Zeeman field B. The local Coulomb interaction 
is U and the interaction between electrons on the two sites is 
U' . The system is coupled to noninteracting leads by hopping 
amplitudes t L / R . 



the approximate fRG treatment of the Coulomb interac- 
tion specific to the present model are discussed. For later 
reference we also give a brief account of the appearance 
of the Kondo ridge for a single dot within our approxi- 
mation scheme. In Sec. IIVI we discuss our results. We 
first describe the effects of the SOI on the spectrum of 
the noninteracting isolated double dot. Next we present 
our results for the linear conductance G(V g ,B) consider- 
ing the entire parameter space. Our work is summarized 
in Sec. 



II. MODEL 

Our multi-level quantum dot model is realized by a 
serial double dot each having a single spin-degenerate 
level (at B = 0) as sketched in Fig. [TJ The Hamiltonian 
of the isolated dot contains several terms 
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with dj a being the creation operator of an electron on 
the dot site j = 1,2 (Wannier states) of spin a, contains 
the conventional hopping t > 0, and the on-site energies 



-l/a 



V g ±5 



(3) 



which can be tuned by an external gate voltage V g . The 
difference of the on-site energies is parametrized by the 
level splitting 5. The effect of SOI is taken into account 
by an imaginary hopping amplitude of spin-dependent 
sign. It is the lattice realization of a Rashba SOI result- 
ing from spatial confinement The Rashba hopping 
term with amplitude a > reads 



-ffsoi = 



T > di. 



H.c. 



(4) 



with the third Pauli matrix a z . Choosing the z-direction 
in spin space for the direction of the SOI breaks the spin- 
rotational invariance. The Zeeman field can have a com- 
ponent parallel to the SOI (that is in z-direction) and one 
perpendicular to it. For this we here choose (without loss 
of generality) the x-direction such that the (local) Zee- 
man term reads 



H z = E [djA a ^w' d i,°' sh 

a,a' .3=1,2 

+ dA a {(T x ) C j. < j>d h < jl cos9 . 



(5) 



For 9 = ±7r/2 the SOI and the -B-field are (anti-)parallel. 
In this case the conventional hopping and the SOI can be 
combined to an effective hopping and the SOI does not 
have a specific effect. In particular, it does not destroy 
the finite- -B Kondo ridges (see Sec. IIVB The local 
Coulomb interaction is included by 
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for the local U > and nearest-neighbor U' > interac- 
tions respectively, with n Ji(T = cr 4j> an< ^ n j = S CT n j,<r- 
By subtracting 1/2 from nj a in the definition of H lnt the 
point V g — corresponds to half-filling of the double dot 
even in the presence of Coulomb repulsion. 

The dot Hamiltonian is supplemented by a term de- 
scribing two semi-infinite noninteracting leads, which we 
here model as one-dimensional tight-binding chains 



lead 



■ E EE 

f}=L,Rj=Q a 



[ C l,j+h 



H.c. 



(7) 



with lead operators c^ +1 a and equal band width 4t. In 
the following we choose r as the unit of energy: r = 1. 
The dot-lead couplings are given by a tunnel Hamiltonian 

#coup = (*Z-4,a C i,l,ff + *H c H,l !< r d 2,<' + H - C J ( 8 ) 
a 

with tunnel barriers set by tL/R. We here consider the 
so-called wide band limit (see e.g. Ref. [l9l ) in which the 
tunnel barriers only enter in combination with the local 
lead density of states (at lattice site 1) evaluated at the 
chemical potential by Tl/r = Trt L ^ R p\ eai is- 



III. METHOD: FUNCTIONAL RG 

A. Flow equation for the self-energy 

We briefly review the applied approximation scheme 
which is based on the fRG method^ focusing on aspects 
specific to the present model. A detailed description of 
the implementation for quantum-dot systems in the ab- 
sence of SOI is provided in Ref. Recent extensions 
including SOI involve inhomogeneous quantum wires ^ 
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Starting point of the fRG scheme is the bare (U = U' = 
0) Matsubara frequency propagator Go of the double dot. 
The leads are projected onto the dot sites and enter via 



the hybridization Tl/ r ^ In the basis 

{|l,t>,|U>,|2,t>,|24}} 
of single-particle states its inverse reads 



(9) 
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with T L / R (<jj) — T L / R sgn(w). This propagator is re- 
placed by one in which low-energy degrees of freedom 
below a cutoff A are suppressed 

Go(iu) = Q(\u\ - A)g (iw) . (11) 

The cutoff A is later sent from oo down to 0. Inserting 
Gq in the generating functional of the one-particle irre- 
ducible vertex functions, an infinite hierarchy of coupled 
differential equations is obtained by differentiating the 
generating functional with respect to A and expanding it 
in powers of the external fields. Practical implementa- 
tions require a truncation of the flow equation hierarchy. 

We restrict the present analysis to the first order in the 
hierarchy and only consider the flow of the one-particle 
vertex, that is the self-energy. It was previously dis- 
cussed analyticall y 19 ' 21 that for a single dot the resulting 
flow equations capture the appearance of a Kondo ridge 
in G(V g ) of width ~ 1.5 U. The accuracy can be im- 
proved by including the flow of the static part of the 
two-particle vertex (effective interaction) which gives a 
Kondo plateau of width ~ U in good agreement with the 
exact result . 19 ' 21 Here we are not interested in such quan- 
titative improvements and instead keep the bare vertex. 
This has the advantage that the resulting flow equations 
for the matrix elements of the self-energy have a sim- 
ple structure which not only allows for a fast numerical 
solution but provides the opportunity to gain analytical 
insights. 

The flow equation for the self-energy reads^ 

A E A j0= l £ ^e- 0+ ^ b ,(^)r Q ', fc ';a, h , (12) 
' ^ uj=±A b,b' 

where the indices a,a',b,b' label the quantum numbers 
(j, a), T a i^ :a j, is the two-particle vertex, and the in- 
teracting Green function G is determined by the Dyson 
equation 

G A {iu)= [GoHiu)-^]' 1 ■ (13) 

The initial condition for Ao — » oo is S A ° = 0. 19 In the 
lowest-order scheme the two-particle vertex T a > b> ia ,b is 
given by the bare anti-symmetrized interaction. As the 



bare vertex is frequency independent, the approximate 
self-energy turns out to be frequency independent. Dy- 
namical contributions are generated only at higher or- 
ders. As the latter are important for the conductance at 
temperatures T > the current approximation scheme 
is restricted to T = 0. The correct temperature depen- 
dence of the (single-dot) Kondo ridge is only captured if 
the flow of a frequency dependent two-particle vertex — 
leading to a flowing frequency dependent self-energy — 
is kept . 22 ' 23 Within our approximation, the matrix ele- 
ments of the self-energy at the end of the flow E A=0 can 
be interpreted as interaction-induced renormalizations to 
the noninteracting model parameters such as the SOI and 
conventional hopping amplitudes, as well as the on-site 
energies and the magnetic fieldJ^ 

B. Computation of the linear conductance 

From the self-energy obtained at the end of the flow at 
A = 0, the full propagator including interaction effects 
is determined via the Dyson equation (|13|) . From this 
various observables can be computed J£ Here we concen- 
trate on the linear conductance. At T = current-vertex 
corrections vanish and the Kubo formula for the spin- 
resolved conductance assumes a generalized Landauer- 
Biittiker form2i 

G a , a . - ^|7;^(0)| 2 , (14) 
h 

with the effective transmission T a ,a> (0) evaluated at the 
chemical potential. For the present setup the transmis- 
sion is given by the (1, <r; 2, a') matrix element of the full 
propagator leading toi£ 

G a ^ = e -r^ L T R \G XM ^t ■ (15) 
h 

C. Single-level quantum dot 

Before analyzing the serial double dot, for later ref- 
erence we briefly discuss the single-level quantum dot 
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within our approximation scheme i 19 i 21 In this case the 
flow equation for the effective (flowing) level position 
V A = V a + oB + is 



dA 



V£ 



(A 



(V, A ) 2 



(16) 



at V g = 0. The renormalized magnetic field B rcn at A = 
hence shows the same exponential behavior (with pref- 
actor B instead of V g ) as the renormalized level position 
Eq. (| 1 T[) except for the reversed sign in the exponent. 
This leads to a dramatic increase of the renormalized 
field B lcn compared to the bare one. For U > the total 
conductance at V„ = 



Vn + aB, a 



and 



with initial condition V A ~°° 

r = r L + r R . 

At B — the level position is spin independent 
V A = V A and the differential equation can be solved 
analyticallyj 19 ' 21 For U 3> T (in the Kondo regime) and 
\V g \ < 0.77 U the solution at A = is 



G(V g = 0,B) = 



2e 2 4Tr.r 
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in the B direction becomes exponentially (set by the 
Kondo scaled ' 19 ' 21 ) sharp, instead of Lorentzian-like of 



It is this exponential pinning of the renormalized level 
position to zero (the chemical potential) which leads to 
the Kondo plateau in the total conductance G(V g , B) — width r for U — 0. In fact, this holds not only at V g = 
G^(V g ,B) + Gi(V g ,B) given by but for all gate voltages within the B = conductance 

plateau. Qualitatively the Kondo ridge of a single-level 
dot in the V g -B plane looks similar to the surrounding of 
one of the B = Kondo ridges appearing in our double 
dot model as e.g. shown in Fig. [2] 



G(V g ,B = 0) = — 



2e 2 4r L r i 



r 2 v 2 + r 2 



(18) 



For U = the level position V — V g is unrenor- 
malized and G reaches its maximum value G max = 
(2e 2 /h) 4r L r^/r 2 only at the resonance voltage V g = 
(Lorentzian resonance of width T). For U > 0, V is 
pinned to zero around V g — for a gate- voltage range of 
width ~ 1.5 U, where the Kondo plateau develops in the 
conductance G ~ G max - 

We next consider a finite magnetic field JS / 0. The 
flow equation for the effective B A — (V A — V A )/2 reduces 
to 



dA (A 



UB a /tt 



r) 2 + {B A y 



(19) 



IV. RESULTS 



A. The noninteracting isolated double dot 



For a detailed understanding of the Kondo ridges it 
is instructive to first study the noninteracting isolated 
double dot. In the basis of Eq. (|9]) the single-particle 
Hamiltonian h^ot is represented as a complex 4x4 matrix 



/ Vg + S + B sin 6 



hdot 
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-t — ia 
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—t + ia 


V g - S + B sin t 
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—t — ia 
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V g -5-Bsint 



(21) 



The corresponding eigenvalues are 

A = V g ± V ' B 2 +t 2 + S 2 + a 2 ± 2B\Jt 2 + 6 2 + a 2 sm 2 
= V g ± JJb~± ^ +S 2 + a 2 sin 2 ^ 2 + a^^cos 2 ^ . 

(22) 

The spectrum is invariant under the transformation B — > 
—B and symmetric at V g — 0. A finite on-site energy 
S > yields an effective hopping t e g = \Jt 2 + 6 2 . For 



vanishing SOI (a = 0) Eq. (|2"2"|) reduces to 

X = V g ± (tea ± B) . (23) 

Naturally, the ^-dependence drops out. For 9 = ±n/2 the 
SOI and the Zeeman field are (anti-)parallel and a can be 
absorbed into an effective hopping t c g = Vt 2 + S 2 + a 2 . 
The eigenvalues are then of the a = 0-form Eq. (|2"3")> . 

For the appearance of a spin Kondo effect (after turn- 
ing on T L / R as well as U and U') it is necessary that 
two degenerate levels of opposite spin are located at 
zero energy (the chemical potential). Zero energy spin- 



degenerate levels are obtained at B = and V g — ±i e ff 
(bonding and anti-bonding states). This will lead to the 
standard spin Kondo effect related to the presence of a 
(spin) Kramers doublet when U and U' are switched on. 
Increasing B the spin-up level of the bonding state and 
the spin-down level of the anti-bonding one approach 
each other. For either a — or 9 — ±ir/2 they be- 
come degenerate (cross) at zero energy for B c — ±i c ff 
and V g = 0. Besides the two B = Kondo ridges devel- 
oping for all dot parameters, for a = or 9 = 7r /2 one 
might thus expect the appearance of two finite- B Kondo 
ridges. Viewed on the basis of the many-body energies of 
the isolated interacting double dot this situation corre- 
sponds to the crossing of a two-particle S = 1, S z = — 1 
state with a corresponding S = 1, S z = state. The 
resulting Kondo effect can thus also be referred to as a 
singlet-triplet one. For a > and an arbitrary angle 9 
between the SOI and the Zeeman field the finite- B level 
crossings turn into avoided crossings and no zero-energy 
degeneracies are possible at finite B. 



B. Numerical results for the conductance 

We next present results for the linear conductance 
G(V g ,B) obtained by numerically solving the flow equa- 
tion (fT2"|) . For our purposes it is sufficient to consider 
only a single set of U, U' and t. The results depend only 
quantitatively on the strength of the local Coulomb in- 
teraction and the inter-dot hopping t, as long as U/T L / R 
and t/tL/R are sufficiently large. We here focus on 
U = U' = 1 andi = l. 



1. Vanishing SOI 

We start our discussion with the case of vanishing SOI. 
This is trivially realized for a = 0. As discussed above, 
for 9 = tt/2 the SOI can be absorbed into an effective 
hopping. In this case the Zeeman field and the SOI are 
(anti-)parallel and only a single Pauli matrix enters the 
Hamiltonian. The physics is thus similar to the a = 
case. 

In Fig. [5] we show G(V g , B) for a left-right symmetric 
setup with t L i R — 0.3 and 8 = in the absence of SOI 
(a = 0). The conductance shows the expected two pairs 
of B = and V g — Kondo ridges. Due to the renormal- 
ization of the inter-dot hopping the B = plateaus are 
not centered around V g = ±i e ff = t = 1 but renormalized 
gate voltages. For the present parameters the renormal- 
ization of the inter-dot hopping and the magnetic field 
(almost) cancel each other such that the finite B ridges 
are located at B c s» ±t e g = t = 1. For U' < U, B c 
is renormalized to smaller values. All Kondo plateaus 
exhibit the same maximal height of 2e 2 /h. This can 
be understood from transforming the dot-lead coupling 
Eq. © into the basis of bonding (b) and anti-bonding (b) 
states of the noninteracting isolated dot. The couplings 
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FIG. 2. (Color online) Conductance G(V g ,B) for t = 1, U = 
U' — 1, t L / R = 0.3, and <5 = in the absence of SOI (a = 0). 



between the leads and the two states are 



b,L/R 



Vt 2 + S 2 T 8 
2y/t 2 + 8 2 



L/R 



(24) 



and T h L/R 



Tj, R/f.t 2 /n/tr 



b,R/L l l/r/ l r/l- ^ 0r ^ ne considered 
left-right symmetric case and (5 = the bonding and 
anti-bonding states have the same total coupling r b /j = 
+ ^b/b ,r. an d the couplings are left-right symmet- 
ric implying unitary conductance. Increasing \B\ from 
B = the Kondo ridges are suppressed on the exponen- 
tial Kondo scale and resonance peaks of height e 2 /h and 
width r develop. The position of the resonance peaks 
varies linearly with V g . Eventually the peaks correspond- 
ing to the spin- up level of the bonding state and the spin- 
down level of the anti-bonding one merge with the finitc- 
B Kondo ridges. Figure [3] shows G(V g ) at fixed different 
B > to further exemplify this. The shoulders appearing 
in the B = bonding (anti-bonding) state Kondo ridges 
are linked to the presence of the anti-bonding (bonding) 
one. Increasing the Coulomb interaction has the two ob- 
vious effects of broadening the Kondo ridges and increas- 
ing the distance between the centers of the plateaus; the 
latter also holds for increasing t. In the lower panel of 
Fig. |3] the respective dot fillings are shown for various 
values of B. For B = the dot occupation exhibits 
plateaus at odd fillings. Their width corresponds to the 
Kondo ridges observed in the conductance. Similarly, for 
a finite magnetic field B — 1 the plateau around V g = 
in the conductance is reflected in the filling. 

We note that G(V g ,B) is symmetric with respect to 
B —B (see the eigenvalues Eq. and V g -> —V g . 

While the former symmetry is given by the Hamiltonian 
and holds for all parameter sets, the latter is specific to 
the parameters chosen here (t^ — £r and 8 = 0). The 
conductance remains symmetric under V g — > —V g if at 
least either 8 = or ti, = t R holds. In these cases the four 
Kondo ridges do no longer reach the unitary value 2e 2 /h 
but exhibit an equally reduced conductance plateau as 
the bonding and anti-bonding states are coupled with 
the same asymmetry and the same total coupling to the 
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FIG. 3. (Color online) Gate-voltage dependence of the con- 
ductance G and dot occupation (ni + n%) for different values 
of B and the same parameters as in Fig. 
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FIG. 4. (Color online) Conductance G(V 9 ,B) for t = 1, {/ = 
{/' = 1, a = with asymmetric couplings to the leads fz, = 
0.3, tfi = 0.5 and finite level splitting 8 — —0.3. 



2. Effect of the SOI 



leads. 

The most interesting situation arises if the left-right 
symmetry and the bonding-anti-bonding state symme- 
tries are simultaneously broken. This is achieved for 
ti, 7^ tit and (5^0. In this case the V g — > —V g sym- 
metry of G(V g , B) is broken as shown in Fig. [U For our 
parameters the couplings of the anti-bonding state have 
a strong left-right asymmetry leading to a B = Kondo 
plateau with significantly reduced conductance (around 
V g = 2.2). The bonding state has a weaker asymme- 
try such that the B = Kondo ridge centered around 
V g = —2.2 almost reaches the unitary conductance. The 
total coupling is larger than Tb- With the breaking 
of the Vg — > — V g symmetry, manifest already from the 
comparison of the two B = Kondo ridges, the finitc- 
B Kondo ridges (centered around V g = 0) are no longer 
necessarily parallel to the V^-axis. In fact they are bent 
with respect to this axis as becomes apparent from Fig. 2] 
It turns out that the direction of bending is always away 
from the state with stronger total level-lead coupling. For 
our model this is the state with larger left-right asymme- 
try. This bending cannot be predicted considering the 
isolated double-dot even by including the Coulomb in- 
teraction as it results from a level renormalization asso- 
ciated with the dot-lead couplings. Related level renor- 
malizations are discussed in Refs. and H^. In exper- 
iments on multi-level quantum dots left-right symmetry 
is difficult to realize and the states at different energies 
will have different level-lead couplings. Therefore finitc- 
B spin Kondo ridges appearing in measurements are ex- 
pected to be generically bent with respect to the B = 
ones. This result is relevant for the understanding of 
finite- B Kondo ridges observed in conductance measure- 
ments on multi-level carbon nanotube quantum dotsi^l 



For the discussion of the effect of the SOI on the Kondo 
ridges we focus on the symmetric case with — t# and 
6 = 0. Figure [5] shows G(V g ,B) for a = 0.6 and an 
angle 8 = 1.56 close to the parallel configuration with 
9 = ir/2. Thus the effective SOI given by the compo- 
nent perpendicular to the direction of the Zeeman field 
is small. For the noninteracting single-particle levels of 
the isolated double dot this implies a small minimal dis- 
tance between the levels avoiding the crossing at finite B. 
Therefore remnants of the finite- B Kondo ridges are still 
observable. The inset shows G at V g — as a function of 
and B. Increasing the SOI component perpendicular 
to the direction of the Zeeman field by deviating from 
6 = 7r/2 obviously destroys the finite- B Kondo effect as 
expected. This provides a way to probe the presence 
of SOI in multi-level dots: after observing finite- B spin 
Kondo ridges one can probe their robustness by chang- 
ing the direction of the magnetic field. Remarkably, the 
SOI introduces an intriguing angular dependence even 
for a magnetic field coupling exclusively to the spin de- 
gree of freedom.— From the angular dependence of the 
gap opening in the single-particle energy spectrum the 
strength of the SOI can be extracted by spectroscopic 
measurements. 



C. Analytical insights 

We next provide analytical insights to our findings in 
the absence of SOI by analyzing the fRG flow equation 
(fT2"]) . To simplify the analysis we focus on the case of 
purely local Coulomb interactions with U' = 0. Due 
to the absence of Fock terms t remains unrenormalized, 
but this only yields quantitative changes compared to the 
results shown in the last subsection. 

In the absence of SOI 8 does not play any role and 
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G/(e*/h) 



can be chosen arbitrarily. For = ±ir/2 the matrix (fT0|) 
is block diagonal and we can restirct the analysis to a 
single spin sector. Introducing V A a = Cj + aB - 
full propagator including the self-energy reads 

<vM = 



£A the 



(25) 



1 



-iT RS ga{uj)-V 2 A 



-t 



-iT LS gn(uj) -V A 



D A (iuj) V ~* ^ 

with the determinant 
D A (iu) = 

[iLo + iT R sgn{io) -V£][iu+iT LSm {w) -V&] - t 2 



(26) 



FIG. 5. (Color online) Conductance G(V 9 ,B) for t = 1, U = 
U' = 1, t L/R = 0.3, and S = with 6 = 1.56 and finite SOI 
a — 0.6. Inset: Conductance G(6,B) at V g —0. 



The zeros of L> A=0 (0) for T L/jR = determine the zero- 
energy levels. For degenerate levels these are responsible 
for the development of the Kondo ridges in the conduc- 
tance. 

Inserting the above expression in Eq. (fT2"|) . the flow 
equations for the local potential are 



U 



dk Vl ° 7T [V A V£ 



At/A 



v 2 A -t 2 )v A 



(A + T R ) 2 V A 



(A + T L )(A + T R ) - *a]a + [(A + T R )V A + (A + I^Ap 



(27) 



with an analog equation for V A a with (1 O 2) and (i f-s- 
i?). 

We first analyze the symmetric case for 5 = (ei = 
£2) and = r# = r. The condition for zero-energy 
levels V 2 — t 2 — implies degenerate solutions for either 
B = and V = ±t or B ren = ±i and V = 0, where the 
renormalized level position or Zeeman field replace the 
bare ones with respect to the noninter acting case. The 
flow equation for the local potential (|27|) reduces to 



d v a _ UV^ 

dA ° IT 

uv A 



{v- A f — t 2 + (A + ry 



[(t/A) 2 -t 2 -(A + r) 2 ] 2 + 4(i/A) 2 (A + r) 2 



-Re 



(A + r + it) 2 + (V- A ) 2 



(28) 



I 

splitting of the central plateau for V in two plateaus of 
half width for V ranging from V g = ±t to ±(t + 0.77 U). 
These are shifted to larger values in Fig. [3] as a conse- 
quence of the renormalization of i. The same substitution 
around the B = Kondo ridges V* = tsgn(V A ) + <tB a 
leads to an exponential suppression of the renormalized 
field as described by Eq. (|19p for the single-level dot. A 
similar behavior in the V g -B plane is found for the finitc- 
B Kondo ridges. Here the renormalization of B leads to 
a shift of the position. Due to the enhancement of the 
renormalized magnetic field, the finite- B Kondo ridge de- 
velops at a reduced field for which B rcn = ±t. This effect 
is superposed by the renormalization of t in Fig. [3] and 
hardly visible. 



resembling Eq. (|16|) for the single-level quantum dot, ex- 
cept for the presence of the finite inter-dot hopping t and 
a factor 1/2 in the definition of T. At B = the flow 
equation for V A — V A — tsgn(V A ) is characterized by an 
exponential pinning to zero for \V\ < 0.77U, inducing a 

I 



We now consider the general asymmetric situation. 
With Vjy 2< 7 = i $ an d introducing the asymmetry 
parameter /3 for the coupling to the leads T L / R — T ± /3, 
the flow equation for V A is 



d v*= u 



V A [{V A f -6 2 -t 2 + f3 2 ] + (A + T)[(A + r)^ A - 206} 



dA 



[(V A ) 2 - 5 2 - t 2 + /3 2 - (A + T) 2 ] 2 + 4[(A + r)vA - ps\* 



(29) 



Analogously, a flow equation for the renormalization of following as it turns out not to significantly affect the 
8 can be derived, which will not be considered in the results. For S ^ and symmetric couplings to the reser- 



8 



voirs /3 = (Tl = Tr), the above equation reduces 
to (f2"5)l . Kondo plateaus of width 0.77 U are hence ob- 
tained in correspondence of degenerate energy levels at 
V = iteff f° r B — and at Vg = for B Ten = ±i e ff- For 
asymmetric couplings to the leads (0 ^ 0), the above Eq. 
(|29p in proximity of the degenerate energy levels can be 
simplified to 



U 



% A - 206 /V 



—V A = 

dA a 7T (A + r) 2 + 4(v£ - ps/vy 



(30) 



in the limit of small asymmetric couplings <C T. The 
terms proportional to /35/T induce a shift in the effec- 



tive level position responsible for the V g 



-V„ symme- 



try breaking. For V g > an effective magnetic field of 
strength B c + 206 /T is required to compensate for the 
shift (corresponding to a larger effective field for 06 > 
in Fig. [3J . This implies that the finite- B Kondo plateaus 
in the V g -B plane are not parallel to the B — axis for 
06 + 0, see Fig. gj 



V. CONCLUSION 

We studied the linear conductance of a serial quan- 
tum dot as a minimal model to describe the influence of 
the SOI on the spin Kondo effect in presence of a Zee- 
man field. 28 Without SOI, the linear conductance as a 
function of an applied gate voltage exhibits characteris- 
tic Kondo plateaus at finite magnetic field B, in addition 
to the ones at B = 0. Interestingly the finite- B Kondo 
ridges are bent with respect to the V^,-axis if the left- 
right symmetry and the symmetry between the bonding 
and anti-bonding states are broken simultaneously. This 
finding is of importance for the understanding of mea- 
surements of the linear conductance of multi-level quan- 
tum dots as a function of Vg and B. In the presence 
of SOI the finite- B Kondo ridges disappear; in contrast 
to the ridges at B = 0, they are not protected by time- 
reversal symmetry. 
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